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Abstract. In this article we study the possibility that, due to non-linear 
couplings, unstable g-modes associated to convective motions excite stable 
oscillating g-modes. This problem is of particular interest, since gravitational 
waves emitted by a newly born proto-neutron star pulsating in its stable g- 
modcs would be in the bandwidth of VIRGO and LIGO. Our results indicate 
that nonlinear saturation of unstable modes occurs at relatively low amplitudes, 
and therefore, even if there exists a coupling between stable and unstable modes, 
it does not seem to be sufficiently effective to explain, alone, the excitation of the 
oscillating g-modes found in hydrodynamical simulations. 



PACS numbers: 04.30.Db,04.40.Dg,97.10.Sj, 97.60. Jd 

1. Introduction 

It is well known that proto-neutron stars (PNS) born in the aftermath of a 
gravitational core collapse are, during the first tens of seconds, convectively unstable 
[1] . In addition, dynamical simulations of core collapse have shown that a considerable 
fraction of the energy emitted in gravitational waves (GW) is actually due to this 
convective motion, rather than to the violent bounce Recently, a new mechanism 

for core collapse supernova explosion has been proposed, based on the excitation and 
sonic damping of core g-mode oscillations an d the same mechanism has also 

been pointed out to be relevant for GW emission [7]. The process begins being 
an advective-acoustic oscillation, and in [5] the authors suggest that the primary 
agent of the explosion is the acoustic power generated in the inner turbulent region; 
the accreting PNS is a self-excited oscillator, and the core oscillation acts as a 
transducer converting accretion energy into acoustic energy. Surprisingly, and despite 
the complexity of the non-linear hydrodynamical effects involved in these simulations, 
the growth times, oscillation frequencies, and temporal evolution of the modes are 
qualitatively and quantitatively similar to the predictions of a linear analysis done 
in [S]; there, the frequencies and damping times of g-modes were computed for an 
evolving PNS, modeled as sequences of quasi-stationary configurations. 

These results put the g-modes in a somewhat different perspective with respect to 
the preceeding literature on the subject; indeed, previous studies of g-mode oscillations 
in neutron stars had shown that their contribution to the gravitational wave signal 
was likely to be too weak to be detectable by current or near- future detectors [5]-[T5]. 

Motivated by these considerations, in this paper we study the spectrum of stable 
and unstable g-modes of a newly-born PNS, and we investigate how effective the 
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coupling is, due to the non linearity of the equations of hydrodynamics, between 
unstable g-modes and oscillating g-modes. The aim is to understand whether 
this coupling induces g-mode oscillations which may significantly contribute to 
gravitational wave emission during the early phases of a proto-neutron star's life. 
To this purpose, we integrate the equations of relativistic stellar perturbations for 
the model of PNS considered in [8], which describes the quasi-stationary evolution 
of a PNS during the first minute after core collapse; we look for the eigenmodes the 
frequency of which, at some point of the evolution, becomes imaginary or acquires 
a negative imaginary part, a behaviour which indicates a mode instability. We then 
discuss to what extent unstable and stable g-modes are coupled. 

It should be mentioned that the existence of an unstable branch of g-modes is 
a natural consequence of the PNS thermal and chemical profiles, since some regions 
of the star are convectively unstable while other regions are not. Thus, even at the 
linear level one can see that both stable and unstable g-modes are present. The linear 
analysis allows to determine which is the fastest growing mode and the frequencies of 
the oscillatory modes. 

The PNS we consider, at very early stage has a radius of about 100 km, and it 
gradually contracts in a few seconds. It has a fixed baryonic mass of 1.6 M Q , which 
corresponds to a final gravitational mass (after deleptonization and cooling) of 1.46 
Mq. The radius of the PNS at t=50 ms after bounce is 50 km but just 1 second 
later, due to neutrino losses from the envelope, the star contracts to a radius of less 
than 20 km. The final radius of the cold, catalyzed NS (after 30-40 seconds) is 12 km. 
In computing the g-mode eigenfrequencies at selected values of time, we have taken 
the corresponding radial profiles (mass, pressure, temperature, lepton content, etc.) 
provided by the evolutionary models of [1], already used in [5]. 

Evaluating the effect of nonlinear couplings between different modes is a 
formidable task; however an estimate of the leading order contributions can be 
provided using a Newtonian approach based on the amplitude equation formalism 
developed in [14] and [IB] . 

The paper is organized as follows. In Section [5] we describe the mathematical and 
numerical setup of our calculations; in Section Owe discuss the structure of stable and 
unstable g-modes in PNS; in Section [4] we study the coupling between different modes 
in terms of appropriate integrals of the radial eigenfunctions; concluding remarks are 
in Section El 

2. Equations and numerical setup 

The equations governing the perturbations of spherical stars in general relativity have 
been derived within different approaches by many authors [E]-[22]. Here we use 
the Lindblom-Splinter (LS) equations derived in [21], which describe dipolar (1 = 1) 
perturbations, and the Lindblom-Detweiler (LD) equations derived in 20J, for I > 2 
perturbations. The metric and the fluid four- velocity of the spherical background 
describing the unperturbed star are 



where the functions v(r), A(r) satisfy the TOV equations (see for instance 23J). We 
shall assume that matter in the star is a perfect fluid. 




- e v dt 2 + e x dr 2 + r 2 (d9 2 + sin 2 0d(f> 2 ) 
e- y/2 , 0,0,0), 



(1) 
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The metric perturbations and the Eulerian fluid perturbations are expanded in 
spherical harmonics Y lm {9, <p), and Fourier-expanded as described in HIJ- In 
particular, the radial component of the Lagrangian displacement is 

Sr = Sr lm (cj,r)Y lm (9,(f>)e iuJt . (2) 

The perturbed Einstein equations 

5Gfj, v = SirST^ 

provide a system of linear, ordinary differential equations (ODE) for the radial 
components of the harmonic expansion of the metric and fluid perturbations. Since 
the structure of the equations for I = 1 and for I > 1 is quite different, we will treat 
the two cases separately. 

I = 1 perturbations are not associated to gravitational wave emission. 
Therefore, since we are neglecting further dissipative effects like viscosity or heat 
transport, oscillation modes are undamped, normal modes and the corresponding 
eigenfrequencies, uik, are real. Dipolar oscillations are described by the LS equations, 
which form a linear, third order ODE system [2"T] . In order to find the mode 
eigenfrequencies, these equations have to be integrated from the center (where 
regularity conditions must be imposed) up to the stellar surface, for different values 
of uj. The solutions satisfy the required boundary condition at the surface, i.e. the 
vanishing of the Lagrangian perturbation of the pressure, only for a discrete set of 
values of w. these are the normal mode eigenfrequencies. The equations depend on 
uj quadratically, and since all terms in these equations and all boundary conditions 
are real, there is no restriction in assuming that the perturbation functions and uj 2 
are real, i.e. in looking for normal modes. Unstable modes correspond to solutions 
belonging to uj 2 < 0. They are non oscillating, exponentially growing modes, with 
growth time r given by 

- = V /= ^- (3) 

T 

For / > 2, stellar perturbations are associated to gravitational wave emission, 
therefore the corresponding modes, said quasi-normal modes, are damped and the 
eigenfrequencies are complex 

uj = 2tiv + - ■ (4) 

T 

v is the mode frequency, and r is its gravitational damping time. These perturbations 
are described by the LD equations, which form a linear, fourth order ODE system |20j. 
These equations are integrated from the center, imposing regularity conditions, to the 
stellar surface, for different values of complex uj. For each u> there exist two regular, 
independent solutions which have to be matched on the stellar surface to allow the 
vanishing of the Lagrangian perturbation of the pressure. This constraint reduces the 
number of independent solutions to one, modulo an overall scaling factor. 

Outside the star, the variables associated to the fluid motion vanish and the 
equations reduce to a second-order ODE for the Zerilli function Z lm (uj,r*) [24] 

^+^Z lm (oj,r*) = V l (r)Z lm (uj,r*) , (5) 

where V (r) is the Zerilli potential and r* = r + 2Mlog(^g — l) is the tortoise 
coordinate. Far away from the star, the Zerilli function behaves as a superposition of 
ingoing and ingoing waves: 

Z lm {r)^Z^^' +Z^ t e-^ . (6) 
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Quasi-normal modes are the solutions of the LD+Zerilli equations for which the 
ingoing component of the asymptotic solution ©, Z^, vanishes. The discrete set 
of frequencies corresponding to these solutions are the mode eigenfrequencies. 

Like the LS equations, the LD equations have real coefficients; however the 
outgoing wave boundary condition Z\™ — cannot be satisfied by a real Z lm . 
Therefore, in order to represent a physical mode the harmonic components of the 
perturbations must belong to complex frequencies. For I > 2 unstable modes are 
solution of the perturbed equations with < 0; an unstable mode corresponds 

to a perturbation which, while oscillating, is increasingly growing in amplitude; the 
oscillation frequency v and the growth time r(> 0) of such unstable mode are given 
by w = 2t:v — i/r. 

The relevant equations for I > 1 have been integrated using an adaptive stepsize 
fourh-order Runge-Kutta method. For I > 2 the solution has been continued outside 
the star by integrating Zerilli's equation §5§ using the continued fraction method [25] ; 
the complex eigenfrequencies have been found as zeroes of the function Z l ™(io), using 
a Newton-Raphson method. 

3. Normal and quasi-normal modes of proto-neutron stars 

Before discussing our results, we shall briefly summarize what is known in the literature 
about the mode structure of non rotating stars in Newtonian gravity. This will serve 
as a guideline to understand our results for relativistic PNS. 

The normal mode spectrum of non-rotating stars has been studied in great detail 
in Newtonian gravity (see for instance |26j.|27j). The main results can be summarized 
as follows. 

• Since Newtonian stars do not emit gravitational waves, in absence of viscous 
forces the square of the mode eigenfrequencies, v\, are real for any value of I. 
The index k denotes the order of the mode. If i/f > 0, the mode corresponds to 
a stable oscillation, while if vf, < the mode is a pure exponential growth. 

• Mode eigenfrequencies group in distinct classes, identified on the basis of the 
restoring force which dominates in bringing back to the equilibrium position the 
generic fluid element displaced by the perturbation; they are named t^-modes if 
the restoring force is buoyancy, p^-modes if it is due to pressure gradients, g- 
mode eigenfrequencies are smaller than those of p-modes, and the two classes are 
separated by the single fundamental mode (/-mode) frequency. 

• When 1=1, the /-mode frequency is identically zero: the /-mode does exist only 
for I > 2. I — 1 g- andp-modes have an interesting characteristic: they correspond 
to a displacement of the geometrical center of the system, i.e. Sr(0) / 0; as 
pointed out in [5H], in the I = 1 case the geometrical center does not coincide 
with the center of mass of the star, which is not displaced. 

• There can exist two classes of g-modes: the gj^ modes, which are stable, and 
the gi \ which are unstable. Whether both branches of modes do exist or not 
depends on the sign of the Schwarzschild discriminant 



— If A = throughout the star (as it is for barotropic EOS), then there are no 
g-modes (or, more precisely, they are all degenerate to zero frequency). 



.4 




d Pj s P 



(7) 
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— If A < throughout the star (convective stability), there are only <7W-m0d.es. 

— If A > throughout the star (convective instability), there are only g-- 
modes. 

— If A < in some region of the star, and A > in some other region, then 
both </ + )-modes and </ - )-modes are present. 

Thus, the sign of the Schwarzschild discriminant reveals the presence of convective 
instabilities, and of unstable g-modes. p-modes and the /-mode are always stable. 

• Typically, for - and j»fc-modes the number of nodes in the eigenfunction of 
(6r)k is equal to k (not counting the node at r — for I ^ 1). The /-mode 
eigenfunction has no nodes inside the star. However, such simple prescriptions 
may not hold in some cases [2"B] . 

Let us now consider our relativistic PNS. In order to find the frequencies of the 
quasi-normal modes, we have integrated the equations that describe the perturbations 
of a spherical star taking as a background the "quasi-stationary" evolutionary models 
derived in pQ. "Quasi-stationary" means that the stellar interior is described by a 
sequence of equilibrium configurations, which have been shown to adequately describe 
the evolution of a PNS for t > 0.1 — 0.2s after bounce [HES]. These models have been 
used in [5] to compute how the quasi-normal mode frequencies of stable, quadrupole 
modes evolve in time. Here we consider the equation of state named GM3, with no 
quark matter in the core, and explore the time interval from t — 0.2s up to t = 25s 
after bounce. We find both stable and unstable g-modes. Unstable modes are present 
because the relativistic Schwarzschild discriminant [30] 

S M = §-(!),§ < 8 » 

is negative in some region inside the star, revealing the presence of convective 
instability^. 

To compare our results with the literature on Newtonian stars (in particular, with 
Figure 1 of [26]), we define a real, square frequency, v 2 , as follows. 

• For I = 1, u% is always real; it is positive for stable modes, negative for unstable 
modes. Therefore we define vfc = (wfc/27r) 2 . 

• For I > 2, Wfc is complex. For stable fluid modes, in general \$loJk\ ^> \^su>k\, 
whereas for unstable modes \^5oJk\ \$tuk\- Therefore, we can define 

(2-rrUk) 2 = (Kcjfe) 2 ~ uj 2 for stable modes 

(2-Ki/k) 2 = — (Sajfc) 2 ~ ui\ for unstable modes . (9) 

The mode structure of the considered stellar model is shown in Figure [1] where v 2 
is plotted for different values of the harmonic index v% is plotted for the first two 
unstable g-modes, g[~^ and g^~\ respectively, for the first two stable g-modes, gj + ' 
and #2 , an d f° r the fundamental mode. The figure refers to t — 2 s after bounce. 
We see that v 2 has a dependence on I similar to that shown in the Newtonian case 
[2"6] . We remind that the main difference is that for I > 1 in the relativistic case io\ is 
complex: the frequency of stable modes acquires a small imaginary part corresponding 
to the damping induced by gravitational wave emission, while that of unstable modes 
acquires a small real part. 

I Notice that the sign convention for the relativistic discriminant JHJ is the opposite to that assumed 
for the Newtonian discriminant 0. 
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Figure 1. Mode structure of a PNS at t = 2s after bounce. 

Let us now focus on the behaviour of unstable g-modes. The lowest dipole mode, 
(i.e. g[ ' for I = 1), at t = 2 s after bounce has a growth time r (-) = 0.51 ms. This 
is very interesting because, since r (_> is extremely small, the amplitude of this mode 

9i 

may have enough time to grow significantly before being damped by other dissipative 
effects, or by non-linear couplings. For I > 2, the growth times of the lowest unstable 
g-modes are even smaller; for instance, for I = 2, r (-) = 0.3 ms. In Figure [2] we show 
how the growth time of the lowest unstable g-modes for I — 1,2 changes during the 
first 25 s of the PNS life. Wc sec that, after an initial decrease, all r <-> 's increase with 

9\ 

time; however, they remain quite small, being their maximum values of the order of 
a few ms. After about 25 s, the star becomes convectively stable and unstable modes 
disappear. 

For I > 2 unstable g-modes have a small oscillating part, which is of the order of 
v ~ 10~ 8 Hz and can be neglected. 

4. Mode couplings 

The existence of rapidly growing, unstable g-modes indicates, at a linear level, the 
onset of an effective convective instability; this motivates our further investigation on 
mode coupling, to see whether unstable modes, coupling to stable ones, may trigger 
gravitational wave emission at frequencies which fall in the bandwidth of ground based 
interferometers, as shown in [5]. 

A first thing one can do is to see whether the corresponding cigenfunctions overlap: 
if they do not, we would have a strong indication that the coupling is ineffective. In 
Figure [3] we compare, as an example, the radial cigenfunctions of the lowest, stable 
and unstable g-modes belonging to I — 1,2. We see that there is indeed an over lap! 
however, this information is not sufficient to infer how effective the transfer of energy 
from convective motion to oscillatory modes can be. 

§ Notice that the eigenfunctions of the stable modes are defined mostly where S(r) > 0, while the 
eigenfunctions of the unstable modes are defined mostly where S(r) < 0. 
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Figure 2. Growth time of the lowest unstable g-modes during the first 25 seconds 
of the PNS life. 




star 



Figure 3. Eigcnfunctions of the lowest stable and unstable g-modes, plotted for 
I = 1,2 as functions of the normalized radius. The Schwarzschild discriminant 
S(r) is positive on the left of the dashed vertical line, negative on the right. This 
picture refers to t = 0.5s after bounce. 



In order to quantify the effectiveness of the couplings, we have used the amplitude 
equations formalism described in [14] and |15j . where couplings between different 
modes have been analysed in the framework of Newtonian gravity. In this approach, 
the Lagrangian displacement Sr(t,r) is expressed as a superposition of modes 

<fr(i,r) = ^ a i (<K W (f), (10) 

i 

where (r) are the mode eigenf unctions, normalized by imposing £jT« (R s tar) / Rstar = 
1, and a % {t) are the corresponding amplitudes. For non-rcsonant modes [jj, the non- 
11 Resonant modes are those for which the condition . kjVj ~ is satisfied, with kj small integers 
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linear coupling between different modes is described by the equation 



Ql 



J\3 



,aV) 2 , 



(11) 



where the frequency lui is in general complex, i.e. u>i — 2^Vi ± — , if some dissipative 
process is in action. In our case, dissipation is due to gravitational waves and, as 
mentioned before, for stable modes we can assume 2irv ^> - while for unstable modes 

The evaluation of the coupling coefficients Qhj , the expressions of which are given 
in |14j and [15] . is a very hard task. Since we need an order of magnitude estimate, 
we shall compute only a leading term. In particular, we shall use the expressions of 
Qljj given in [15 , but we have checked that using those given in [14] the results are 
of the same order of magnitude. The leading term we consider is 

3(2-^)4,. 



Ql j:j - iC, 



(12) 



where 



Z\jj = I dQY^YjYj 



A = drpr 2 \& 



C] 



(13) 



p j V p 

Yi is the spherical harmonic corresponding to the i-th mode, and p,p, Sp, Sp are 
the density, the pressure, and their Lagrangian perturbations. If all quantities are 
expressed in geometrical units (G = c = 1) as powers of km, then Q\jj has dimension 
km -1 , which can be converted to s _1 multiplying by c. 

Let us first compute the coefficients Q\ u , i.e. the self-coupling coefficients (see 



eq. pip ), for an unstable mode. Since in this case 



Qm 



2 m ji 



eq. ([T2]) gives 



1, 2, we find 



(14) 



and, for the lowest unstable g-modes corresponding to I 

\Q l=1 \ -3.0-10 13 s- 1 

\Q l=2 \ - 1.3 -10 7 b _1 . 

Assuming that Q\ u < 0, i.e. that the effect of self-coupling is that of saturating 
the mode, using the values of and |Q /=2 | found above, we can estimate the 

saturation amplitude as the value of a 1 for which a 1 

1 



0, i.e. 



(15) 



We find: 



8.2 cm 
150 m . 



for 1 = 1 Sr(R star ) 
for I = 2 5r(R star ) 

Next we want to check whether unstable modes can trigger stable modes, eventually 
making them unstable. Therefore, we need to evaluate the coefficients which describe 
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Table 1. Dimensionlcss coupling coefficients QL»(i| at ) between stable modes 
(index i) fi and unstable modes (index j) 
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4- icr 2 



the coupling between the i-th stable mode and the j'-th unstable mode. If the z-th 
mode is stable, eq. (fTT|) becomes 

a 1 ~ i2^a* - -a 1 + Ql vl (a z ) 3 + Q\ jja \a j f . (16) 

The purely imaginary terms on the right-hand side of (116|) determine a shift of the 
frequency Vi and are not relevant for the mode stability; therefore, we are interested 
only in the real part of the right-hand side of (fTB)) . If we define 

571 



I lien eq. (fTTj) can be written as 



[l + Q^V) 2 ] . (18) 

We must then evaluate the dimensionless quantity ;(a J sat ) 2 , where a 3 sat is the 
saturation amplitude of the unstable mode j, previously computed. In particular, 
we want to check whether the sign of the term in square brackets in equation (fT%|) can 
be reversed. 

It should be noted that the stable i-th mode must belong to / > 2, because 
modes belonging to I = 1 are not associated to gravitational wave emission, and for 
them 1/Tj = 0. This means that the mechanism of mode enhancing we are studying 
is triggered by gravitational wave emission, as in the case of the well-known CFS 
instability [3T] , 

In Tabled] we show the values of \Q\jj | (a-*) 2 ; the index j refers to unstable modes 
g^, \ where n = 1, 2 is the order of the mode, and I = 1, 2; the index i refers to stable 
modes gty and /; (fundamental mode), with 1 — 2. In all cases we find 

IQy(0 2 «i; (is) 

this means that the term in square brackets cannot change sign and that the effect 
of unstable modes on stable modes is too weak to induce a significant growth of their 
oscillation amplitudes. 

5. Concluding remarks 



We have determined the frequencies of stable and unstable modes of a newly born 
proto-neutron star during the early phases of its life. Our goal was to understand 
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whether the coupling between stable and unstable modes could provide a mechanism 
to convert the energy of convective motion into (/-mode oscillations, enhancing 
gravitational wave emission. 

Our results indicate that the existing coupling do not seem large enough to 
explain, alone, the excitation of the g-mode oscillations found in hydrodynamical 
simulations. Thus, other processes should be considered as driving the excitation of 
these modes, as, for instance, accretion from the mantle [4], [32] . 
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